Critical exponents of the driven elastic string in a disordered medium 
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We analyze the harmonic elastic string driven through a continuous random potential above the 
depinning threshold. The velocity exponent (3 = 0.33(2) is calculated. We observe a crossover 
in the roughness exponent £ from the critical value 1.26 to the asymptotic (large force) value of 
0.5. We calculate directly the velocity correlation function and the corresponding correlation length 
exponent v = 1.29(5), which obeys the scaling relation v = 1/(2 — £), and agrees with v FS , the 
finite-size-scaling exponent of fluctuations in the critical force. The velocity correlation function is 
non-universal at short distances. 



Driven elastic manifolds in disordered media model the 
physics of systems as diverse as charge density waves Q , 
interfaces in disordered magnets 0, contact lines of liq- 
uid menisci on rough substrates |3[, vortices in type-II 
superconductors Q and crack propagation in solids H . 

An elastic manifold, driven through disorder by an ex- 
ternal force, undergoes a dynamic phase transition 
that arises from competition between the driving force 
and the pinning energy due to the disorder, mediated by 
the elasticity of the manifold. Analogous to an equilib- 
rium phase transition, the driving force acts as the con- 
trol parameter, and the average center-of-mass velocity 
v of the manifold acts as order parameter. Two phases 
of respectively zero and non-zero order parameter are 
separated by the critical force / c . At forces below this 
depinning threshold / c , the disorder 'pins' the manifold 
and for long enough times, the velocity of the manifold 
is zero. Above / c , the manifold continues to advance 
in avalanches. When the threshold is approached from 
above (/ — + f£), the mean velocity tends to zero, and 
typical length, width and duration of the avalanches di- 
verge. This critical divergence is characterized by two 
independent scaling exponents. 

Much effort 0, IE IS E3| has been spent on calculat- 
ing the universal exponents, particularly for driven elas- 
tic manifolds in the limit of quasi-static motion. In this 
limit, inertial terms are neglected and the force acting on 
the manifold is assumed independent of velocity. The net 
force on the manifold comprises a constant driving force 
/, a position-dependent random force r\ and an elastic 
restoring force. Choosing a harmonic short range elastic 
force leads to the following equation of motion for the 
manifold h at zero temperature 

d t h(x, t) = f + r)[x, h{x, t))] + dlh{x, t). (1) 

In general, the manifold is represented as a single- valued 
function h(x, t) defined over a D-dimensional transversal 
space x, moving in a D + 1 dimensional disorder. In this 
paper, we consider one discrete dimension D = 1 with 
periodic boundary conditions for the string of length L 
(x + L is identified with x) and for the disorder of lateral 
extension M (h + M is identified with h, up to a winding 
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FIG. 1: Mean velocity v of the elastic string vs / — / c smpl , 
(L = 128, M = 195, 18 samples). The window of critical 
power law behavior is sandwiched between the single-particle 
regime v ~ f 1/2 at very small /-/ c smpl — fitted to eqU — 
and linear v ~ / behavior at large /. Inset: critical window 
for larger system sizes L = 512, . . . , 2048 with M ~ -L . Slope 
on log- log plot yields the velocity exponent j3 = 0.33(2). 



term) . The continuous disorder of unit strength and unit 
correlation range is constructed like in |ll| . 

In what follows, we numerically investigate the dynam- 
ics of the 1 + 1-dimensional string above the depinning 
threshold. Convergence of the dynamical solution is en- 
sured by exploiting the particular analytical structure 
[13, Il3| of the equation of motion: it has a unique pe- 
riodic solution for each disorder sample in the t — > 00 
limit. Having found the periodic solution, we are thus 
certain to have reached the asymptotic regime and to 
have shaken off all influence of arbitrary initial condi- 
tions. The asymptotic periodic solution is constructed 
to desired precision using a continuous integration rou- 
tine for eq. ^ Averaging observables over one period and 
over disorder samples, we calculate the velocity expo- 
nent (3 and describe in detail how the correlation length 
diverges as /— >/+. 

In the thermodynamic limit, the string velocity obeys 
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a power law: v ~ (/ — / c )' 3 for / — > / c . On finite sys- 
tems, the critical force /^mpi fl UCTlU ates from sample to 
sample, putting a limit on how small the control parame- 
ter /— (/c mpl ) can be made without introducing undesired 
corrections to scaling relations; a limit from which previ- 
ous numerical calculations suffer. 

We are able to determine the exact sample-dependent 
depinning threshold f* mp thanks to a recent algorithm 
j 1 11] . The sample critical force itself shows non-negligible 
fluctuations of the order of <jf c « L~ 1 / l/FS . The algo- 
rithm also finds the final critical configuration h c and 
the roughness exponent ( at depinning pL3J | - 

Knowing the critical force /™pi Q £ sample, we 
plot the time- and disorder-averaged velocity v against 
/ — /c" 1 • Thus we eliminate the statistical noise due to 
the fluctuating critical force, and obtain extraordinarily 
clean data. Furthermore, the control parameter f—f^ mpl 
can be made arbitrarily small, which is not possible when 

using/- (/ c smpl ). 

The mean velocity on a finite sample shows three dif- 
ferent regimes (see figure^: For very small / — /,? mp \ the 
motion of the entire string is correlated, and it behaves 
effectively like a single particle (v ~ (/ — J™ 1131 ) 1 / 2 ). At 
intermediate forces, the string is correlated on length- 
scales £ and moves collectively (v ~ (/ — /c mp1 )^)- At 
high forces the motion is essentially uncorrelated (v ~ /). 

At / = /™ pl , the finite dynamical system (eq. ^1 (in 
the t — > oo limit) undergoes a saddle node bifurcation 
between a static (pinned) and a periodic (depinned) so- 
lution; the global velocity minimum changes from sta- 
ble to unstable. For very small positive / — f^ mpl , the 
string spends the major part of its time-period passing 
through the velocity minimum, and negligibly little time 
completing its orbit. The motion through the minimum 
is dominated by the mode h with the largest eigenvalue 
in a linear stability analysis around the critical config- 
uration h c . The mode h moves at a velocity equal to 
/ — /® mpl plus quadratic and higher corrections in h. 
Hence we can express the time spent inside the mini- 
mum as: / f = / f _ f ^ l+ch2 = H(/ - /r p T 1/2 , to 
first order in / — /™ lpl . The remaining time T spent 
outside the velocity minimum depends only weakly on 
/ - / c smpl . This yields for the velocity as a function of 
force: 

M/v =T = T + Ti(f- y^pi)-Va. ( 2 ) 

As shown in figure^ this function — combining the effec- 
tive single-particle exponent one half (v ~ (/ — /£ nip1 ) 1 ' 2 ) 
for very small forces with the saturation at slightly larger 
forces — fits the data perfectly. 

At intermediate forces, the dynamics leave the single- 
particle regime and enter the regime of critical collective 
motion. The cross-over takes place when the dynamic 
correlation length £, which diverges as £ ~ (/ — f c )~", 
equals the system size L 0]. ft follows that the dif- 
ference between the corresponding cross-over force and 



the sample critical force scales as L^ 1 / 1 ' . This scaling 
is quite slow because of the small exponent 1/v w 0.7, 
as our data confirms. The other cross-over, between the 
critical region and the linear regime, is independent of 
system size. 

Between the two cross-overs lies the window of collec- 
tive behavior, small and slowly growing with L: Only 
strings of length L > 512 allow to see any significant ev- 
idence for critical behavior, and even for L = 2048 (see 
inset of figure Q the window is less than two decades. 
These pronounced finite-sample-size effects limiting the 
critical window (also observed in the CDW model flq]') 
could obscure critical behavior in experimental situa- 
tions. For example, in an experiment of a liquid-solid 
contact line advancing on a disordered substrate 0] the 
sample size, i.e. the capillary length acting as upper cut- 
off, is L w 200 in units of the disorder correlation length. 
For these small samples we expect the window of collec- 
tive behavior to be hardly noticablc and finite-size effects 
to dominate. 

When extrapolating data to the thermodynamic limit, 
we have to carefully choose the lateral sample size: M 
has to be of the order of the typical width W of a string 
of length L at depinning. When increasing L, W scales 
as L^, £ being the roughness exponent. Hence the lat- 
eral sample size M has to scale as Z>, too. Otherwise, 
if M is scaled with an exponent £' < £, the periodic 
sample is too short, the string wraps around it, and gen- 
erates correlations which mix in properties of the charge 
density wave model (CDW . M ~ 1), itself member of a 
different universality class . If on the other hand M 
scales with £' > £, the sample is too long and contains 
~ M/W ~ independent critical configurations of 

size L x W . Each of these configurations has a slightly 
different local critical force, and the critical force of the 
entire sample is given by their maximum, and not their 
mean. Consequently the critical force of a sample of size 
M W overestimates / c , even in the thermodynamic 
limit. 

We therefore use a scaling of M ~ L/> , with strings of 
length L = 512, 2048, when analyzing the mean ve- 
locity. Inside the window of critical collective behavior, 
we find (3 = 0.33(2) (see inset of figure This value is 
consistent with 0.2(2) from a two-loop functional renor- 
malization and with 0.35(5) from a simulation of an 
advancing magnetic domain wall in a 2D random Ising 
system It is larger than previous estimates on con- 
tinuous systems 0.22(2) 0| and on automaton models 
0.25(3) In these studies, however, v was plotted 

against / — (/c mpl ), and j3 was probably measured partly 
inside the critical window and partly inside the finite-size- 
effect region where the velocity saturates. This naturally 
leads to f3 being underestimated. 

At the depinning threshold the typical size £ of 
avalanches diverges. Whereas previous work did not suc- 
ceed in accessing £ directly, we are able to compute this 
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FIG. 2: Connected velocity- velocity correlations C v (x) for 
different values of f-fc, (L = 512, M = 20000, 700 samples). 
Inset: correlation length £ as obtained from a fit to C v (x) = 

dynamic correlation length from the velocity correlation 
function and, independently, the structure factor. In our 
continuous calculation, contrary to automaton models, 
the instantaneous velocities v(x,t) are readily available, 
and the velocity correlation function can be computed 
unequivocally: 

L 

C v (x) = {v(x)v(0)) c = J2(( v ( x '+ x > t )-v){v(x',t)-v)) t . 

x' 

Two points on the string taking part in the 
same avalanche have correlated instantaneous velocities. 
Therefore, the typical length £ of avalanches shows up as 
the characteristic length scale in the velocity correlation 
function. 

Analyzing this correlation function C v (x) we find a di- 
verging correlation length £ (see figure|2J). In order to ex- 
tract £, we perform for every given / — f c a two-parameter 
(£,k) fit of the correlation function to the functional form 

C v (x)~C Q -x- K e- x/ t. (3) 

The inset of figure |21 shows the correlation length £ ob- 
tained as a function of f—f c . It diverges as £ ~ {f — fc)~ u 
with v = 1.29(5). This value agrees with fpg — 1-33(1) 
found in a finite-size-scaling analysis of sample-to-sample 
fluctuations of the critical force plj . Hence we confirm 
that the dynamic correlation length exponent and the 
fmite-size-scaling correlation length exponent are identi- 
cal v = i/ FS Furthermore, v o bey s the statistical 
tilt symmetry relation v = 1/(2 — £) 10] . 

The velocity correlation function is found to be non- 
universal: the exponent n characterizing the algebraic 
decay at small distances depends on / — f c , unlike previ- 
ously thought 0, Q ■ This is why the correlation function 
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FIG. 3: The structure factor (S(q)} ~ l/q 1+2C for different 
values of f - fc (L = 512, M = 20000, 200 samples). Two 
regimes of different roughness £ are visible: At small q-values 
C, = 1/2, while for large q, £ = 1.26. The crossover tends 
toward zero for f f c , illustrating the diverging correlation 
length. Curves for different / — f c are offset along the y-axis, 
lines are guides to the eye. 

C v (x) for different values of / — f c cannot be plotted on 
one single master curve as a function of x/£. We have 
not attempted to establish a definite functional relation 
between k and / — f c . 

The non-universal behavior of the velocity correlation 
function at small distances might be linked to the fact 
that in one dimension the harmonic model's roughness 
exponent is greater than one: The average distance be- 
tween neighboring points grows without bound, imply- 
ing that any real string described by this model would 
eventually break |22j] . Preliminar y da ta on a model with 
higher than harmonic elasticity [23j does indicate uni- 
versal velocity correlations, suggesting that the harmonic 
short range model is indeed exceptional. 

An avalanche of typical length £ has a typical width 
w which scales like w ~ £*». As a second method to esti- 
mate the roughness exponent £, we study the time- and 
disorder-averaged structure factor S(q), which behaves 
like l/q 1+2 '', when defined as S(q) = (h(q)h(—q)) c with 
%)=£*e^(x). 

At large driving forces / — > oo, the roughness expo- 
nent of the harmonic elastic string is C = 1/2: The noise 
due to the disorder becomes equivalent to thermal noise, 
which can be seen from expanding the disorder term r\ 
in eq. ^ of motion - transformed to the center of mass 
reference frame - in powers off -1 « yielding a 8{t) 
correlated noise to first order. At depinning, the rough- 
ness exponent takes on the value £ — 1.26(1) 

Inside the critical window, both these values for £ show 
up in the structure factor, and characterize two different 
ranges of q separated by the inverse of the correlation 
length: for large wave- vectors (q > =?) the structure 
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factor displays the critical roughness, whereas at small 
wave- vectors (q < ^p) it shows the thermal roughness. 
Figure El tracks the crossover between the two regimes as 
/ — » /+, illustrating the diverging correlation length £. 
The corresponding correlation length exponent v is con- 
sistent with the value calculated from the velocity corre- 
lation function, but less precise. 

In our calculation of the correlation function, the lat- 
eral size M must be chosen sufficiently large for the string 
not to notice the periodicity of the disorder — other- 
wise, the structure factor displays at small q, in addi- 
tion to the two regimes mentioned, the roughness of a 
CDW (Ccdw = 3/2 [l3|). This means that the small- 
est g-modes have to completely decorrelate on the scale 
M. The mode q decorrelates in a time ~ q~ z with z the 
dynamic exponent, z not only controls the typical du- 
ration of the avalanches (of size £(/ — f^Y), but also the 
decorrelation time of all length scales |24j . in particular 
the string length L. Note that inside the critical window 
L > £. This implies that at a mean velocity of order 
unity, i.e. at the upper end of the critical window, the 
smallest mode needs a distance of roughly L z to decorre- 
late. But z w 1.5 is larger than £. We therefore have to 
choose an M ^> Z> in our analysis of S(q) and C v (x), in 
order to eliminate auto-correlation effects in the largest 
length-scales. 

On the contrary, in our initial analysis of the mean 
velocity (see figure 0), we were justified in adopting a 
scaling M ~ Lt because the effective string viscosity — 
which determines /3 — is caused by pinning on length 
scales below £. The periodicity M ~ L 1 * thus does not 
influence the velocity exponent p. 

At very large M, the sample critical force f^P 1 be- 
comes too large to be used as the critical force. We there- 
fore calculate a mean critical force from a finite-size- 
scaling ansatz: The average sample critical force is taken 
to depend on the sample length L as / c °° - (/ c smpl (L)) ~ 
£-1/Vfs. The asymptotic value / c °° = 1.913(2) is then 
used as a mean critical force when analyzing the velocity 
correlations and the structure factor. 

When investigating the correlation length £ directly 
by means of C v {x) and S{q), we do not need to know 
ysmpi j n or d er ^0 avoid the cross-over from the critical 
window into the finite-size-effect region: We simply limit 
our analysis to sufficiently small values of £/L. In addi- 
tion, the smallest value of / — f£° analyzed is larger than 
typical fluctuations in f^ mpl , which are therefore negligi- 
ble. 

In conclusion, we analyze the quasi-static dynamics of 
the harmonic elastic string driven through a random po- 
tential, above the depinning transition. We calculate the 
exact periodic solution for each sample, and encounter 
large finite-sample-size effects, which will have repercus- 
sions in experimental situations. Knowing exactly the 



sample critical force enables us to identify the limited 
window of critical collective behavior and to determine 
the velocity exponent (3 = 0.33(2). We investigate the 
velocity correlation function, determine the correlation 
length exponent v = 1.29(5), and confirm that it both 
obeys the statistical tilt symmetry and agrees with the 
finite-size-scaling correlation length exponent. Surpris- 
ingly, we find a non-universal functional form for the ve- 
locity correlation function: the exponent k describing the 
short distance behavior depends on the control parame- 
ter / - f c . 
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